Time evolution of dynamic propensity in a model glass former. 
The interplay between structure and dynamics. 
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By means of the isoconfigurational method we calculate the change in the propensity for motion 
that the structure of a glass-forming system experiences during its relaxation dynamics. The re- 
laxation of such a system has been demonstrated to evolve by means of rapid crossings between 
metabasins of its potential energy surface (a metabasin being a group of mutually similar, closely 
related structures which differ markedly from other metabasins), as collectively relaxing units (d- 
clusters) take place. We now show that the spatial distribution of propensity in the system does not 
change significantly until one of these d-clusters takes place. However, the occurrence of a d-cluster 
clearly de-correlates the propensity of the particles, thus ending up with the dynamical influence of 
the structural features proper of the local metabasin. We also show an important match between 
particles that participate in d-clusters and that which show high changes in their propensity. 

PACS numbers: 61.20.Ja; 61.20.Lc 



I. INTRODUCTION 



The determination of the physical basis of the emer- 
gence of glassy relaxation when a liquid is supercooled 
under its melting point represents a subject of intense 
research 1 ^! 3 -!^ 6 -. Striking aspects of such behavior are 
the fact that the dynamic observables change dramati- 
cally with the supercooling (while static quantities show 
at most very mild differences) and that the dynamics 
varies by orders of magnitude from one region of the 
system to anotheri 2 -' 3 -' 4 ^' 6 -. The dynamics of these sys- 
tems slows down very fast in this regime (as temper- 
ature is reduced) and the relaxation has been believed 
to proceed by means of cooperatively relaxing regions 
whose time scales and sizes grow considerably as the 
temperature is decreasecU^i&^&i. This heterogeneous 
scenario has been corroborated by the discovery, both 
experimentally and theoretically, of the existence of dy- 
namical het er ogeneit ies£ & iO) 11 : 12 : 13 : 14 : 15 . Some studies 
identified cooperative motions of a small number of par- 
ticles (a few percent) which move collectively, often in 
a string-like fashion, by a distance close to the particle 
diamete r 11 : 12 ' 13 : 14 . More recently, collective motions of 
a significant fraction of the particles of different regions 
of these systems which form relatively compact clusters 
have been foundi 6 -. These very rapid and sporadic events, 
which were termed "democratic motions" , drive the sys- 
tem from one metabasin (MB) of its potential energy sur- 
face, a group of similar closely related structures 3 - : 16 ' 17 , 
to another while the structural relaxation (the so-called 
a relaxation) is performed by a small number of such 
events. These relatively compact clusters ("democratic" 
clusters or d-clusters^ 6 -) have been identified in differ- 
ent glassy systems like a binary Lennard- Jones system 
and supercooled water— and represent natural candi- 
dates for the cooperatively relaxing regions proposed long 



ago by Adam and Gibbs 7 -. A recent inhomogeneous 
mode-coupling theory of dynamical heterogeneity has re- 
lated them to the (fractal) geometrical structures carry- 
ing the dynamical correlations at timescales commensu- 
rable with that of the a relaxation 1 ^. Additionally, a 
recent experimental and computational work in a glassy 
polymer provided experimental support to the MB-MB 
transitions and d-clusters^ -. 

Even when it would be intuitively expected that the 
dynamical heterogeneities are related to structural het- 
erogeneities in the sample, the determination of the exis- 
tence of a causal link between structure and dynamics re- 
mains as another main open problem in the field 2 - ' 21 : 22 . A 
recent beautiful idea holds the promise to shed some light 
in this regar d 21 ' 22 . By means of the isoconfigurational en- 
semble (IC) it has been determined that the propensity 
of the particles (their tendency to be mobile) is in fact 
determined by the initial configuration of the sample and 
that the particles with the higher propensity are not ho- 
mogeneously distributed within the sample but arranged 
in relatively compact cluster o 21 ' 22 . The time extent of the 
influence of the local structure on the dynamics has also 
been shown to be on the order of the MB residence time, 
a timescale shorter than the a— relaxation time (r a )^. 
In a previous work, it has been shown for the binary 
Lennard-Jones system that the high propensity regions of 
a given initial configuration represent unblocked regions 
wherein d-clusters occur in the subsequent dynamics (for 
any given isoconfigurational realization or trajectory ini- 
tiated in such configuration)^!. 

The aim of the present work is to complete this pic- 
ture by revealing the fact that not only the local structure 
constraints the resulting dynamics but to elucidate the 
role of the d-clusters in reformulating the spatial varia- 
tion of propensity for motion, thus making evident the 
mutual transference of constraints between structure and 
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dynamics. To this end, and by means of molecular dy- 
namics (MD) simulations of the binary Lennard-Jones 
model we shall generate many ICs at different times over 
a given trajectory (thus determining the spatial distri- 
bution of propensity on configurations at different times 
of the trajectory) and we shall show that at the times 
when the d-clusters occur the system shows the greater 
changes in propensity, thus producing a clear propensity 
de-correlation. 



II. MODEL SYSTEM AND METHODS 

A. The binary Lennard-Jones system 

We performed a series of MD NVE simulations for a 
widely used model of fragile glass former: the binary 
Lennard-Jones system consisting of a 3D mixture of 80 
% A and 20 % B particles, the size of the A particles 
being 47 % larger than the B one a 14 ' 16 i 25 . We shall show 
results from systems at temperature T = 0.5, density of 
1.2 and N = 150 particles^. This system size avoids 
the interference of results from many different subsys- 
tems (metabasins) while being free of finite size effects, 
as shown in Ref. [lg. We have also simulated a large sys- 
tem (N — 8000 particles) which yielded an equivalent 
dynamical behavior (however, in order to look for MBs 
and d-clusters, one has to study the behavior of subsys- 
tems of N = 150 and thus, the large system should be 
decomposed in many small subsystems wherein to apply 
the concept of M B 17 ' 26 ). Given this fact, in the rest of 
this work we shall refer to a system of N = 150. At 
low temperatures (close and above the critical tempera- 
ture predicted by the mode-coupling theory of the glass 
transition, T c — 0.435) this system presents dynamical 
het erogeneit ie a 1 ' 1 6 : a small number of particles move 
cooperatively a distance that is comparable to the inter- 
particle distance. These "fast moving" (or "mobile") 
particles are not homogeneously distributed throughout 
the sample but are arranged in clusters usually made 
of string-like groups of particle s 14 ' 16 . The dynamics is 
most heterogeneous at time t* defined by the maximum 

in the non-gaussian parameter a-i (t) , a.2 (t) = ^ \ — 1 , 
which measures the deviation of the self part of the van 
Hove function (the probability at a given time of find- 
ing a particle at distance r from its initial position) 
from a brownian behavior 14 . This quantity is located 
at the end of the [3 - beginning of the a relaxation (the 
crossover from the caging to the diffusive regime in the 
mean-squared displacement (r 2 (t)) plot) and constitutes 
the characteristic time for dynamical heterogeneities (in 
this case t* = 400 for T = 0.50^) . Additionally, t* de- 
pends strongly on temperature and grows quickly as we 
move towards T<*L. However, not all the mobile particles 
within a t* time span contribute decisively to the a relax- 
ation, as we have recently demonstrated^. Instead, the 
a relaxation is driven by a series of a few MB-MB transi- 



tions which are triggered by the occurrence of large com- 
pact clusters of medium-range-mobility particles called 
democratic particlesJ~ Additionally, the mean-residence 
time in an MB has been estimated to be close to t*— . 



B. The distance-matrix method, metabasins and 
d-clusters 

We now describe briefly the distance-matrix method 
to study MB dynamics (see Refs. I161I231 for details): We 
perform an MD simulation and record equally-spaced 
configurations (for example 101 configurations, as in 
Refs . Il6li23h for an a-relaxation total run time (r Q = 
4000, thus consecutively recorded configurations are sep- 
arated in time by <f> = 10 % t* = 40) and build the fol- 
lowing distance matrbsi 6 .^, A 2 (t', t") = ± J2?=i l r i(*') _ 
Ti(t")\ 2 , where rj(t) is the position of particle i at time t 
(since A 2 (£',£") = A 2 (t",t') this is a triangular matrix). 
Thus A 2 (t',t") gives the system- averaged-squared dis- 
placement of a particle in the time interval that starts 
at t' and ends at t" . In other words, this distance 
matrix contains the averaged-squared distances between 
each recorded configuration and all the other ones. For 
this study (as all studies dealing with MB s 16 ' 17 ' 28 ), we 
must investigate small systems, since for large systems 
the results originated from different subsystems would 
obscure the conclusion o 16 i 17 ' 28 . Thus, we used 150 parti- 
cles. However, we also found the same qualitative results 
for small subsystems immersed in a big one, thus rul- 
ing out the possibility for finite size effects (we repeated 
the study for subsystems of 150 particles within a large 
system of 8000 particles, that is, we focused only on a 
portion of the large system) . 




t' /100 



FIG. 1: Contour plot of the distance matrix A 2 (t',t") for 
T — 0.5. The gray levels correspond to the values that are 
given to the right. From Ref. [l^ . 

Figure [1] shows a typical behavior for trajectories with 
T = 0.5. The gray level of the squares in A 2 (t',t") de- 
picts the distance between the corresponding configura- 
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tions, the darker the shading indicating the lower the dis- 
tance between them. Were the dynamics homogeneous 
in time, we would expect a distance matrix with a dark 
main diagonal and a continuous fading as we move away 
from such diagonal. However, from the island structure 
of this matrix a clear MB structure of the landscape is ev- 
ident. That is, islands are made up of closely related con- 
figurations (low A 2 ) which are separated from the con- 
figurations of other islands by large distances. In other 
words, the dynamics of the system is inhomogeneous in 
time. We can estimate the typical residence time in the 
MBs for this T (from island sizes) as qualitatively on the 
order of t* . Given the small system size we expect this 
to be a good estimate (however, this timescale clearly 
depends on system size, since for a large system differ- 
ent subsystems would be undergoing MB-MB-transition 
events at different times). Thus, MB-MB transitions (the 
crossings from one island to another and which last 10 
% to 20 % of t*) are fast events compared to the times 
for the exploration of the MBs. The study of MB-MB- 
transition events has been done previously^, revealing 
the decisive role of large compact d-clusters of medium- 
range- mobility particles. These clusters are responsible 
for the a relaxation (completed after 5-10 such events) 
and represent potential candidates for the cooperatively 
relaxing regions of Adam & Gibbs^. The compact na- 
ture of these events relevant to the a relaxation has been 
shown to be in accord with the geometrical structure of 
the dynamics correlations at large timescales on the order 
of r a within an inhomogeneous mode-coupling theory (at 
variance from the less dense structures, compatible with 
string-like motions, expected at the shorter timescales of 
/3-relaxation time)^. 

The democratic particles that comprise the d-clusters 
that trigger MB-MB transitions were defined as that 
whose mobility was greater than r t h = 0.3 within the 
time interval </> = 40, and its fraction is represented by 
the function m(t, 4>)r^. Thus, for the system size under 
study we found that on the order of 40 - 60 particles were 
involved in a d-cluster— . 

In Fig. [2] we show for the same trajectory and total run 
time interval the function S 2 (t,9), the system-averaged- 
squared displacement of the particles within a time in- 
terval 9 (black-solid curve). This function is defined as: 

5 2 {t,6) = A 2 {t-9/2,t + 9/2) 
1 N 

= ^£l r *(*-<V 2 )- r ^ + <V 2 )l 2 - 
t=i 

Thus 8 2 (t,9) is A 2 (t',t") measured along the diagonal 
t" = t' + 9 and hence the average of this quantity over 
different start times t gives the usual mean-squared dis- 
placement for time lag 6. For this plot we have chosen 
9 = 160, a value that is significantly smaller than the 
a— relaxation time (r Q = 4000) but still sufficiently larger 
than the time of the microscopic vibrations (=0(1)). A 
comparison of this 5 2 (t, 9) with Fig. [T] shows that S 2 (t, 9) 
is showing pronounced peaks exactly then when the sys- 
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FIG. 2: Solid-black line (right ordinate): Averaged-squared 
displacement 5 2 (t,8) for the trajectory given in Fig. [1] Solid- 
dark-gray line (left ordinate): The function m(t,(f>) which 
gives the fraction of democratic particles, i.e., particles that 
moved more than the threshold value rth — 0.3 in the time 
interval [t, t + <fi]. Dashed- light-gray line (left ordinate): The 
function (t, <j>) which gives the fraction of particles with 
a great propensity-change, i.e., particles that changed their 
propensity more than the threshold value £ = 0.21 in the 
time interval [t,t + </)]. Vertical bars (left ordinate): The func- 
tion M(t,(f>) with gives the fraction of particles that are both 
democratic and have £ > 0.21. The values of 8 and <j> are 
respectively 160 and 40. 

tern leaves an MB. Thus we see that changing the MB 
is indeed associated with a rapid motion as measured in 
S 2 (t, 9). In Fig. [51 we have included the fraction of demo- 
cratic particles m(t, 4>)^- as a function of time (dark-gray- 
solid line). The comparison with 5 2 (t,9) shows that the 
fraction of these particles is indeed large whenever the 
S 2 (t, 9) increases rapidly. This fraction is on the order of 
30 % to 40 % of all the particles^ and, thus, significantly 
larger than one would expect from 47rr 2 G s (r, 4>) if one in- 
tegrates this distribution from rth to infinity (which gives 
-0.16). 

C. The isoconfigurational method. The role of the 
local structure: Dynamic propensity 

To calculate the propensity for motion we use the iso- 
configurational (IC) method introduced in Ref. [U In it 
one performs a series of equal length MD runs (trajecto- 
ries) from the same initial configuration, that is, always 
the same structure (the same particle positions) but each 
trajectory with different initial particle momenta chosen 
at random from the appropriate Maxwell-Boltzmann ve- 
locity distribution (that is, one builds an IC ensemble). 
For any given time, each run or trajectory presents cer- 
tain mobile particles or dynamical heterogeneities. How- 
ever, the mobile particles and corresponding clusters of 
mobile particles differ from run to run since the mobil- 
ity of the particles in a single run is not determined 
by the initial configuration^ 2 .. Propensity of a particle 
for motion in the configuration at time t (its tendency 
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to be mobile at the instantaneous time t in the trajec- 
tory given in Fig. [TJ for a fixed time interval of length 
4> is defined as (Ar?(t))ic = Pi(t), where (• • ■ )ic indi- 
cates an average over the IC generated at time t and 
Arf(t) = |r'i(t, <j>) — Yi(t)\ 2 . Here, r^t) is the position 
of particle % in the configuration at time t of the trajec- 
tory given in Fig. Q] and where one of the ICs is started, 
and r'i(t,4>) is the position of the same particle at the 
end of one of such IC trajectories of length <f>. This 
definition of propensity is a generalization of the origi- 
nal definition 2 - since in this work we are interested in 
studing how propensity evolves with time during a given 
MD trajectory and to relate such propensity changes to 
the dynamical events that occur in such trajectory. At 
low temperatures propensities for any given configuration 
at time t are not uniform throughout the sample and high 
propensity particles are confined to certain (relatively 
compact) region a 22 ' 23 . Thus, while particle mobility is 
not reproducible from run to run, the spatial variation 
in the propensity is completely determined by the ini- 
tial configuration, reflecting the influence of structure on 
dynamic s 22 ' 23 . While the mean value of the propensity 
[{Pit)) = N- 1 J2? =1 Pi(t)\ depends on the length of the 
time interval in which it is calculated ((f>) , the spatial vari- 
ation of propensity does not depend on it for times not 
too small (down to 10 % t* or even less) 21 ' 29 . In other 
words, if one calculates the propensity at <fi — t* or higher 
(of course that for very long times larger than r Q when 
the system becomes diffusive, the propensity distribution 
gets uniform) or if one does it for a timescale of 10 % t* , 
the particles with high and low propensity are the same 
in any case. In this work we chose a time = 40 = 10 % 
t* to calculate propensity since at this short time the IC 
trajectories do not have time to abandon the local MB 
(perform a d-cluster) and we can safely say that we are 
sampling the short-time vibrations (without interference 
of the dynamics at larger times). A propensity much 
greater than the mean value, a great tendency to be mo- 
bile, is thus a clear indication that the particle is not 
"comfortable" in its actual position, that is, it is struc- 
turally unjammed. 



III. RESULTS 

With the tools already described (subsection III C[) and 
in order to elucidate the influence that the dynamical 
events (subsection III B[) have on the local structure of 
the system, we generated a single dynamical trajectory 
(Fig. |T|) and started many ICs over it. That is, we stored 
configurations from an MD run each <p time units (a to- 
tal of 101 configurations, thus the run was for a total 
length of r a = 4000) and generated 101 different ICs of 
500 trajectories from each of them. Thus, we determined 
the propensities for each of the 101 equally spaced con- 
figurations over the given MD run (we mention that we 
have obtained similar values of propensity using 200 tra- 
jectories for each IC instead of 500). The way the spatial 



distribution of propensity (the values of the propensity 
for each of the different particles) changes from one con- 
figuration to another, gives an idea of the time evolution 
of the local structural constraints. To quantify how simi- 
lar or different is the propensity of the different particles 
in configurations at times t' and t", we calculate the fol- 
lowing cross-correlation function R = Y^Li Pi> where: 

Ri = [Xi - {X)} [Y {¥)] ^ 

{EttlXi {xw} 1 ' 2 {EtiW - (Y)] 2 } 1/2 

In it, Xi and Yi are respectively Pi{t') and Pj(i"). Pi(t) 
is the propensity of particle i calculated over the IC gen- 
erated from the configuration at time t, in the trajectory 
seen in Fig. [TJ Besides {X) and (Y) are respectively 
(P(t')) and {P{t")). Then, with all pairs of times t' and 
t" we generate the squared matrix R(t',t"). Figure [3] 
depicts R(t',t"), which shows a structure which matches 
that of the MBs of Fig. [TJ We can see that large values 
of R (dark regions) occur at times when Fig. [TJ shows 
a metabasin structure. This means that the propensity 
of the particles within an MB are positively correlated, 
that is to say, the propensity of each particle is similar in 
all the corresponding ICs. Since R is generally small for 
ICs not belonging to the same MB, the propensity of the 
particles at such times do not present a neat correlation. 
Moreover, the borders of the squared islands of Figs. [JJ 
and [3j show certain small negative values, thus indicat- 
ing a slight anticorrelation in propensity once the system 
abandons the local MB (high-propensity particles change 
to low propensity, or intermediate propensity, and vice 
versa). This fact means that MB-MB transitions and 
d-clusters reformulate the propensity pattern, by rear- 
ranging the regions of high, medium and low propensity. 
This is also consistent with a rapid loss of the memory 
of the local structure (the structural constraints of the 
present MB) once the system escapes from the MB. Wc 
also note that if we calculate propensities at larger times 
(say, for example at <fi = t* instead of <f> = 40), the results 
are consistent, and even more conspicuous. 

We recall that the function R(t',t") of Fig. [3] is built 
upon "structural information" while A 2 (i',t") of Fig. 
Q] relies on dynamical facts. R(t',t") is based on the 
propensity for motion of the particles calculated at very 
short-time intervals (j> = 40 (at times t' and t") and thus 
reflects the cross-correlations between the local structural 
constraints of the configurations at the corresponding 
times. 

In order to better quantify the role of d-clusters and 
MB-MB transitions in changing the propensity map, we 
calculate the distribution of (root-squared) changes in the 
particle propensities, £j = \Pi(t + 4>) — Pi(i)! 1 / 2 . A large 
value of £j implies that particle i has significantly changed 
its tendency to be mobile at the time t + (f> with respect 
to its tendency at time t. The average of the distribution 
function of & for all the 100 time intervals (47r£ 2 G s (£, <ft)), 
as depicted by the solid line in Fig. 2J would represent 
an analogue of the self part of the van Hove function 
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FIG. 3: Cross-correlation function R(t' ,t") between propen- 
sities calculated over ICs generated from configurations at 
times t' and t" that belong to the trajectory given in Fig. Q] 
The propensities are calculated for a <f> — 40 time interval. 
The gray levels correspond to the values that are given to the 
right. 



in the dynamics of the system (the probability that a 
particle after a time <f> be found at distance r from its 
original position). If we restrict this calculation to the 
times when the system is within an MB (and thus to 
the squared islands of Figs. [T] and [3|), we get a similar 
curve (as shown in Fig. 2] in line with circles), since the 
system spends most of the time within MBs and thus this 
behavior dominates the whole distribution. However, if 
instead we calculate £ for the time intervals where an 
MB-MB transition occurs (the times when the system 
shows the crossings between islands in Figs. [1] and [3]), 
the distribution (line with squares) is clearly displaced 
to the right, thus indicating that a great enhancement in 
the changes in particle propensity occurs at such times 
(as can also be learnt from Fig. 0| . 

We also calculated the fraction of particles that dis- 
play large changes in £ and denoted this function as 
m^{t,4>). As a threshold we use £ > 0.21 in order to 
have values comparable with the function m{t, 4>) (that 
is / o > 2 G s (r,0)dr = J^ 2i eG s (C,m = 0.16/(4tt)). In 
Fig. [3] we have included this function (light-gray- dashed 
line). Direct inspection of such plot shows the clear cor- 
relation that m$ (i, (/>) exhibits with S 2 (t, 9) and the func- 
tion rn(t,<f>). The fraction of matches (function M(t,<j>), 
vertical bars) between particles that are both democratic 
and that exhibit a large propensity change (in a <f> time 
interval) is also indicated. Finally, to further show the 
good homology between democratic particles (which dis- 
placements are greater than 0.3 in a time interval) and 
particles that change significantly their propensity (par- 
ticles with £ > 0.21 in a time interval) in a typical 
MB-MB transition, we show in Fig. \5\ the spatial distri- 
bution of both kinds of particles for the MB-MB tran- 
sition that occurs in the time interval [680,720] in the 
trajectory seen in Fig. [T] Like the d-cluster consisting of 
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FIG. 4: Radial probability density 47r£ 2 G s (£, <f>) of finding a 
particle with a (root-squared) propensity change of magnitude 
£ in a (j> = 40 time interval. Solid line: average over the whole 
run. Line with circles: average over time intervals within 4 
selected MBs. Line with squares: average over 4 selected MB- 
MB transitions. Selected times are the same as in^. Because 
of the noise in the data, these functions have been smoothed. 



democratic particles (light- and dark- gray spheres), the 
particles with high £ (black- and dark-gray spheres) are 
arranged in relatively compact clusters. Both clusters 
occupy the same region in space and many particles are 
both democratic and also have £ > 0.21 in the [680,720] 
time interval as can be seen in the figure (represented by 
dark-gray spheres). 

As a final remark, we should like to mention that the 
behavior above expounded is characteristic of dynamical 
heterogeneities and thus emerges within the temperature 
interval where these heterogeneities are present. At tem- 
peratures lower than T = 0.5 the behavior is very similar 
(in fact the results are a bit neater since the dynamical 
heterogeneities are more conspicuous) but the timescales 
grow since t* increases with decreasing temperature (400 
for T = 0.50 and 5000 for T = 0.446 27 ). In Figs. H 
and [8] we show the situation for T = 0.446. We mention 
that the functions in Fig.0] (T — 0.50) are very similar to 
those for T = 0.446, and the threshold r — 0.3 for demo- 
cratic particles in a <f> = 40 time interval for T — 0.50 (10 
% t*(T)) is also reasonable in a <p — 500 time interval for 
T = 0.446 (10 % t*(T)). Thus, we kept using the value 
£ > 0.21 for a criterium of high propensity change in a 
(j) = 500 time interval for T = 0.446. 

On the other hand, at higher temperatures (higher 
than approximately T = 0.6), the system looses the dy- 
namical heterogeneities and becomes homogeneous and 
diffusive. Consequently, the distance matrix A 2 lacks 
the island structure (as the non-gaussian parameter ui 
begins to vanish) and the relaxation becomes diffusive, 
without the presence of d-clusters^. Additionally, as 
showed in Ref. [22], the propensities also get uniform at 
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FIG. 5: Cluster of democratic particles (light- and dark- 
gray spheres) and high-propensity change (black and dark- 
gray spheres) for the MB-MB transition in the time interval 
[680,720] in the trajectory seen in Fig. [1] Dark-gray particles 
are both democratic and have a high-propensity change. 



high temperature. 
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FIG. 7: Solid-black line (right ordinate): Averaged-squared 
displacement S 2 (t,9) for the trajectory given in Fig. [6] Solid- 
dark-gray line (left ordinate): The function m(t,cj>) which 
gives the fraction of democratic particles, i.e., particles with 
displacement r > 0.3 in the time interval [t,t + 4>\. Dashed- 
light-gray line (left ordinate): The function m((t,tf>) which 
gives the fraction of particles that changed their propensity 
more than £ = 0.21 in the time interval [t,t + </>]. Vertical bars 
(left ordinate): The function M{t,<j>) with gives the fraction 
of particles that are both democratic and have £ > 0.21. The 
values of 9 and <f> are respectively 2000 and 500. Propensities 
are calculated for a (j> = 500 time interval over ICs of 200 
trajectories. 
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FIG. 8: Cross-correlation function R(t',t") between propen- 
sities calculated over ICs generated from configurations at 
times t' and t" that belong to the trajectory given in Fig. [6] 
The propensities are calculated for a <j> — 500 time interval. 
The gray levels correspond to the values that are given to the 
right. 



IV. CONCLUSIONS 



FIG. 6: Contour plot of the distance matrix A 2 (t',t") for 101 
equally-spaced (500 time units) configurations of a trajectory 
at T = 0.446. The gray levels correspond to the values that 
are given to the right. 



The isoconfiguratiorial ensemble has provided us with a 
means of determining regions of the sample that present 
different tendency for mobility, thus enabling to mea- 
sure the influence that the local structure exerts on the 
dynamic o 21 ! 22 . The presence of relatively compact re- 
gions with high propensity for motion indicates domains 
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made up of particles with high tendency to be mo- 
bile wherein the events responsible for the a relaxation 
occu r 23 ' 24 ' 29 . The fact that the spatial distribution of the 
very short time propensit y 21 ' 29 is sufficient to signal such 
regions and that some of such high propensity particles 
will take part of a d-cluster in the different IC trajectories 
at a later tinuM^ would mean that these particles are 
not at ease (not blocked) in the local structure and try for 
some time without success to relax their condition until 
eventually, a collective motion of many of them is able 
to perform a large scale relaxation event. While a rig- 
orous complete comparison is not possible at this stage, 
this behavior seems to be in accord with an appealing 
very recently proposed soft mode explanation of glassy 
dynamics 3 "^. In such scenario^, the high propensity re- 
gions present at a given configuration and the irreversible 
reorganization regions that occur later on are causally 
correlated with the localized low-frequency normal modes 
of such configuration and which persist for timescales on 
the order of t*. We note that the black islands of Figs. [3] 
and [5] indicate that the spatial distribution of propen- 
sity (the tendency of the particles to move when evalu- 
ated at very short times) persist for times \t' — 1"\ ~ t*. 
Within such timescales, the system changes many times 
the inherent structure (IS, local minimum or basin of at- 
traction in the potential energy surface), an event that 
entails a small localized particle rearrangement, but has 



not been able to perform a collective long-range rear- 
rangement characteristic of a d-cluster since it has been 
confined within the local metabasin, a collection of struc- 
turally very similar configurations. The d-cluster, which 
brings the system out of such MB, occurs within the high 
propensity regio n 24 ' 29 and has been related to the soft 
modes present in the system at such time o 19 ' 32 ' 33 . To 
summarize, in this work we have shown that not only 
the local structure poses its constraints on the dynam- 
ics of glassy systems, but that the opposite is also valid, 
that is, the dynamical events responsible for the a relax- 
ation clearly modify the propensity pattern. Our results 
demonstrate that when a d-cluster (and the correspond- 
ing MB-MB transition) occurs, the spatial variation of 
propensity is reformulated (high- and low-propensity re- 
gions clearly de-correlate), thus making evident the mu- 
tual interplay between structure and dynamics. 
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